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ABSTRACT 

The extended Kepler mission, K2, is now providing photometry of new fields every three months 
in a search for transiting planets. In a recent study, Foreman-Mackey and collaborators presented 
a list of 36 planet candidates orbiting 31 stars in K2 Campaign 1. In this contribution, we present 
stellar and planetary properties for all systems. We combine ground-based seeing-limited survey data 
and adaptive optics imaging with an automated transit analysis scheme to validate 21 candidates as 
planets, 17 for the first time, and identify 6 candidates as likely false positives. Of particular interest 
is K2-18 (EPIC 201912552), a bright {K =8.9) M2.8 dwarf hosting a 2.23 ± 0.25 i?® planet with 
Teq = 272 ± 15 K and an orbital period of 33 days. We also present two new open-source software 
packages which enable this analysis. The first, isochrones, is a flexible tool for fitting theoretical 
stellar models to observational data to determine stellar properties using a nested sampling scheme to 
capture the multimodal nature of the posterior distributions of the physical parameters of stars that 
may plausibly be evolved. The second is vespa, a new general-purpose procedure to calculate false 
positive probabilities and statistically validate transiting exoplanets. 

Subject headings: catalogs — planetary systems — planets and satellites: detection — stars: funda¬ 
mental parameters 


1. INTRODUCTION 


The Kepler telescope (Borucki et al. 20T0| has led to 
a revolution in stellar and planetary astropnysics, with 
7305 “objects of in terest” and 4173 “plaii et candidates” 
discovered to date ( BorucMetalj2011a|b Batalha et al. 


|2'()y||Mullally 

y of this sampli 


et al.||2015[ Rowe et al. 


2013 Burke et al.l 


20151. The fidelity of this sa mple is high: most of these 


candidates are truly planets (IMorton fc Johnson 2011 
Fressin et al.||2013l |Dfeert et al.||2015|. 


i'he mechani¬ 
cal failure ot two reaction wheels on the spacecraft led 
to a repurposing of the spacecraft into the K2 mission, 
in which the telescope points a t fields near the ecl iptic 
plane for ^ 75 days at a time (Howell et al. 20141. In 
this observing strategy, two axes ot motion ot the space- 
craft are controlled by the two remaining reaction wheels, 
while the roll of the spacecraft is balanced with Solar ra¬ 
diation pressure and quasiperiodic thruster firing. As a 
result, the detector drifts relative to the sky at the rate 
of ^ 1" hr“^, with rapid corrections due to thruster fires 
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approximately once every six hours. Over the full dura¬ 
tion of each campaign, the targets remain near the same 
location on the detector but both the slow drift an d the 
corrections are observable by eye (|Barentsen||2015 1. 

K2 light curves produced with aperture photometry 
contain substantial pointing-induced photometric varia¬ 
tions caused by the star’s apparent motion over a poorly 
defined fiat field. Worse yet, these variations occur on 
timescales similar to transit signals, potentially masking 
the observational signature of a planet passing between 
Kepler and its host star. 

There has been considerable effort to recover these 
planetary signals, and to date six planets hav e been con¬ 
firmed orbiting three stars in the K2 data (Armstrong 


et al. 2015 Crossfield et al. 12015 Vanderburg et al.|2015|. 

What is common to all ot tnese methods are that removal 


of systematics is considered a step to be undertaken be¬ 
fore the search for planets. Under this strategy, it is 
implicitly assumed that the systematics are removed per¬ 
fectly, while retaining all of the astrophysical signal. Of 
course, it is impossible to perfectly separate the astro- 
physical and instrumental signal, and such a technique 
is prone to either over-fitting, in which some of the as¬ 
trophysical signal is also removed, or under-fitting, in 
which some of the instrumental systematics remain. A 
better strategy is to simultaneously fit both the signal 
and the systematics, as is common practice in cosmology 
and, increasingly, in radial velocity searches for plane¬ 
tary systems (e.g. Ferreira fc Jaffe |2000[ Boisse et al. 


2()11[ Haywood et al.||2(Jl^ Grun blatt et al.||2015| .' 


h'oreman-iVIackey et al. simultaneously tit both 

the systematics and potential planetary transit signals 
in a search for transiting planets. They assume that 
the dominant trends in the observed stellar light curves 
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are caused by spacecraft motion and are shared by 
many stars. They then run principal component anal¬ 
ysis (PCA) on all stars to measure the dominant modes, 
modeling each star as a linear combination of 150 of 
these “eigen light curves” and a transit signal. This 
method enables fitting without over-fitting, and also per¬ 
mits marginalization over uncertainties induced by the 
systematic model. Therefore, any uncertainties in the 
systematics can be propagated into uncertainties in de¬ 
tected planet parameters, instead of assuming the sys- 
tematics are understood perfec tly. Using this technique, 
Foreman-Mackey et al. (20151 detect 36 planet candi- 
date s orbiting 31 stars in K2 Camp aign 1 data. 

In Foreman-Mackey et al. (2015), only transit prop- 
erties are provided, not absolute parameters about the 
planet or the star. Additionally, the authors follow the 
convention of the Kepler team to include any transit 
event as a candidate system rather than a false positive if 
a secondary eclipse is not detected: there is no enforced 
upper limit on the allowed planet radius. The authors 
intentionally make no effort to separate true transiting 
planets from astrophysical events that mimic the appear¬ 
ance of transits, such as an eclipsing binary (EB) with 
a high mass ratio, similar to the Kepler team’s list of 
“objects of interest.” 

In this paper, we present stellar and planetary pa¬ 
rameters for each system. We also analyze the false 
positive probability (FPP) of each system using vespa, 
a new pub licly ava i lable, general-purpose implementa¬ 
tion of the Morton (2012|) procedure to calculate FPPs 
for transiting planet^ Through this analysis, as well 
as archival imaging, ground-based seeing-limited survey 
data, and adaptive optics imaging, we are able to con¬ 
firm 21 of these systems as transiting planets at the 99% 
confidence level. Additionally, we identify six systems as 
false positives. 

This paper is organized as follows. In Section 2, we de¬ 
velop stellar properties through photometric and spectro¬ 
scopic data. In Section 3, we combine the derived stellar 
properties with K2 data to infer planet candidate prop¬ 
erties. In Section 4, we combine adaptive optics and ra¬ 
dial velocity observations with both archival and modern 
ground-based, seeing limited survey data and an analysis 
of the transit parameters to calculate FPPs. In Section 
5, we discuss potentially interesting systems, including a 
mini-Neptune orbiting an M dwarf which receives a sim¬ 
ilar insolation to the Earth. In Section 6, we summarize 
and discuss our results. 


2. STELLAR PROPERTIES 
2.1. Photometry 

With the exception of one star in our sample (K2-18), 
we do not have spectroscopic data with which to char¬ 
acterize the stellar properties. Additionally, there are 
no measured parallaxes for any of these stars. Instead, 
we rely on photometry. For each system, we query the 
VizieR data base of astronomical catalogs (Ochsenbein 
et al.||2000|). We record the B, V, g', r', and i' mag¬ 
nitudes and their uncertainties from the AAVSO Pho- 
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2Ul^. We also record the J, H, and K mag¬ 


nitudes and their uncertainties as found in the 2MASS 


All-sky Catalog of Point Sources (Cutri et al. 2003) and 
the W1 — W3 WISE magnit udes and uncertain ties from 
the ALLWise Data Release (Cutri et al. 2013). For all 
except two of our targets, the Wi band is only an up¬ 
per limit, and in the remaining two cases, the photomet¬ 
ric uncertainity in VF4 is at least an order of magnitude 
larger than those in Wl — 1F3, so we do not use WA for 
any system. These data are reported in Table 1, and a 
color-color diagram showing the r — J, J — K colors of 
our candidates is included as Fig. [T] 



J-K 


Fig. 1. — Color-color diagram displaying r — J,J — K photometry 
for targets observed by Kepler during the original mission (black), 
with our K2 Campaign 1 planet candidates overlaid (red). Also 
included is the location of the Sun (yellow) and host stars of pre¬ 
viously confirmed K2 planets (blue). 90% of our candidates have 
photometry consistent with later spectral types than the Sun. 


2.2. Stellar Models 

To convert the observed photometric data into physical 
properties for each star, we used the new publicly avail¬ 
able isochrones Python module pq a general-purpose in¬ 
terpolation tool for the fitting of^ella r models to pho - 
tometric or spectroscopic parameters (Morton 2015al. 
This software does trilinear interpolation in mass-age- 
[Fe/H] space for any given set of model grids, thus being 
able to predict the value for any physical or photometric 
property provided by the models at any values of mass, 
age, and [Fe/H] within the boundaries of the grid. 

This enables a set of observed properties {{xi,ai}), ei¬ 
ther spectroscopic, photometric, or both, to define a like¬ 
lihood function to be sampled: 


1 r/n\ 1 ~ (^)) 

\nC{e) oc -- 2^ ^2 


( 1 ) 


where li{0) is the isochrone model prediction of property 
i at the given parameters 0. If the observed properties in¬ 
clude any apparent magnitudes, then 9 includes distance 
and extinction in addition to mass, age, and [Fe/H] . 

In this work, we us e grids from the Da rtmouth Stellar 
Evolution Database (Dotter et al. 2008) at Solar values 
of [a/Fe] =0.0 and helium abundance Y = 0.2741, which 
come packaged with the isochrones modul e. We then 
infer the stellar parameters using MULTINEST (Feroz et al. 


http: //github. com/timothydmorton/isochrones 
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20091, an implementation of a mnltimodal nested sam¬ 
pling algorithm, for each host star conditioned on the 
observed photometric properties as presented in Table 
1. MULTINEST is designed to sample mnltimodal poste¬ 
riors, where other samplers snch as MCMC algorithms 
often strnggle. Given the mnltimodal nature of our pos¬ 
teriors, this scheme is optimal for capturing parameter 
space on the subgiant branch where these stars could re¬ 
side. We include a prior on stellar metallicity represen¬ 
tative of the observed metallicities of stars within 1 kp c 
of the Sun, following the results of Hayden et al. (20151, 
and a Salpeter-slope prior on mass up to the maximum 
mass available in the model grids of 3.7 Mq. 

During the sampling process, we fit for Galactic ex¬ 
tinction as one of our physical parameters. We include 
the WISE bandpasses by applying the relative extinc¬ 
tion values between SDSS, 2MASS, and WISE calcu¬ 
lated by Davenport et al. (20141. In each step of our 
fitting process, we draw a value for Ay, calculate the 
expected extinction in all b andpasses Ax ass uming the 
Ry = 3.1 reddening law of (Fitzpatrick||1999 1, and then 
measure the likelihood of our model stellar fit to the 
observed apparent magnitudes. We apply a uniform 
prior ranging from zero to a maximum extinction value 
of 0.2 and marginalize over extinction in our final de¬ 
termination of stellar parameters. The N ASA/IPAG 
Extragalact i c Dat abase, which reports the SchlaJI^_fc 
Finkb einer (20111 recalibration of the Schlegel et al. 
extinction map as measured by GUHE/DIKHE 


and IK AS/ISSA, suggests that typical Ay extinction val¬ 
ues to the edge of the galaxy at this high Galactic latitude 
are ^0.1 mag, so our upper limit appears to be justified. 

Such a scheme enables us to infer the statistical uncer¬ 
tainties on the mass, radius, and effective temperature. 
However, we are subject to biases induced by system- 
atics in the models themselves. There is some evidence 
that the Dartmouth models may under-predict radii of 
M dwarfs by ~ 15% when compared to o ther methods 


(Newton et al. 2015 Montet et al. 20151. Such an ef- 
fect may be the result ot the Dartmouth model reliance 
on BT-Settl atmospheres, which are based on incomplete 
molecular line lists and have been shown to predict near- 
IR colors that are too blue (Thompson et al.||2014 ). 

As our stellar results are niiodel-dependent, we caution 
users who intend to use these parameters for other works, 
such as exoplanet population studies. When available, 
stellar parameters inferred through other techniques such 
as asteroseismology or spectroscopy should supersede 
these values. We note the observed photometric param¬ 
eters are consistent with spectroscopically derived pa¬ 
rameters for stars with published spectra, and consistent 
with typical mod el-dependent uncert ainties from photo¬ 
metric data (e.g. Huber et al. 20141. We provide full 
samples of our posteriors on the physical parameters for 
e ach starF^ 


Bastien et al. (2014) use the “granulation flicker” in the 
Kepler light curves to suggest that approximately 50% 
of planet host stars have evolved off the main sequence 
onto the subgiant branch, so that both the host stars and 
their planets are larger than previously reported. Simi¬ 
larly, in K2 Gampaign 1 we may expect to find evolved 
stars in a sample of planet candidates, although we may 

http://www.astro.princeton.edu/~tdm/k2/ 


expect the effect to be lessened due to the high Galac¬ 
tic latitude of Gampaign 1. Indeed, we find this to be 
the case. Two stars, EPIG 201257461 and 201649426 are 
definitively evolved stars, with inferred masses less than 
2 MQbut radii above 8 R©. For approximately one third 
of the others, we find the stellar radius posterior dis¬ 
tribution to be bimodal, with both main sequence and 
subgiant models of the stars being consistent with the 
photometric data. This number is consistent with our 
expectations of the number of subgiant contaminants in 
the Gampaign 1 field (K. Stassun 2015, private commu¬ 
nication) . Future observations to measure the parallaxes 
of these stars, such as with Gaia, will be helpful in dif¬ 
ferentiating between these two models to determine more 
precisely the stellar, and thus the planetary, radii. 

2.3. SuperNova Integral Field Spectrograph (SNIFS) 
and SpeX Spectroscopy 

A near-infrared spectrum of K2-18 was o btained using 
the u pgraded SpeX (uSpeX) spectrograph (|Rayner et al. 
20031 on the NASA Infrared Telescope Facility (IRlF'j 
on 2015 January 29 (UT). SpeX observations were taken 
using the short cross-dispersed mode and the 0.3x15" 
slit, which provides simultaneous coverage from 0.7 to 
2.5/rm at i? ~ 2000. The target was observed at two 
positions along the slit to subsequently subtract the sky 
background. Eight spectra were taken following this pat¬ 
tern, which provided a final signal-to-noise ratio (S/N) 
of > 150 per resolving element. The spectrum was flat 
fielded, extracted, waveleng th calibrated, and st acked 
using the Spextool package (Gushing et al. 2004). An 
AOV-type star was observed immediately alter the tar¬ 
get, which was used to create a telluric correction using 
the xtellcor package (Vacca et al.||2003). 


An optical sp e ctrum was obtaine d using SNIFS (Alder- 
ing et al. 2002 Lantz et al. 2004) on the University of 
Hawaii 2.2m tdescope on the nignt of 2015 January 30. 
SNIFS provides simultaneous coverage from 3200-9700A 
at a resolution of ~ 1000. Final S/N of the spectrum 
was > 100 per resolving element in the red (~ 6000A). 
Details of the SNIFS reduction, including dark, bias, and 
flat-held corrections, cleaning the data of bad pixels and 
cosmic rays, and extra ction of the one-dim ensi onal spec¬ 


trum are described in Bacon et al. (2001) and 


et al. (2006). Flux calibrati on was perfo r med using a 


Aldering 


separate pipeline described in Mann et al. (2015). 

Teff was calculated bv comparing our optical spectra 
with the GFIST suit^/5 of the BT-SETTL version of 


the PHOENIX atmosphere models (Allard et al. 2013|, 
which gave a temperature of 3 503 ± 60 K. I Vlore de tails 
of this procedure ar e given in Mann et al. (|2013b ) and 
Gaidos et al. (2014). This method was used because it 
is known to accurately reproduce empirical T^g values 
from long-baseline optical interferometry [Boyajian et al.| 

(l2M^. 

Metallcity was de termined using the procedures from 
Mann et al. (2013a), in which the authors provide em- 
pirical relations between atomic features and M dwarf 
metallicity, calibrated using wide binaries. We adopted 
the weighted mean of the H— and AT—band calibrations, 
which yielded a metallicity of 0.09±0.09. 

http://phoenix.ens-lyon.fr/Grids/BT-Settl/CIFIST2011/ 











































































4 


We combined the derived Tes and [Fe/H ] values with 


the em pirical Teff-[Fe/H]-i?* relation from |Mann et al. 
(20151 to compute a radius. Accounting for measurement 


and calibration errors in [Fe/H] and we calculated a 
radius O.394±O.O38i?0. We use these parameters instead 
of the derived photometric properties for this target, al¬ 
though we note the two are consistent at the Icr level. 

The full list of stellar parameters adopted in this paper 
is included in Table 2. 

3. PLANET PROPERTIES 


In Foreman-Mackey et al. (20151, only parameters di- 
rectly observable from the Kz light curve itself were re¬ 
ported: the period, time of transit center, and tran¬ 
sit depth. With stellar properties now in hand, we 
can convert these observational results into fundamen¬ 
tal parameters of each planet candidate. For each can¬ 
didate, we fit the light curve using a phys ical transit 


model (Mandel & Agol 2002 


ously with a systematics model s i milar to the one de- 


Kipping 20101 simultane- 


scribed by Foreman-Mackey e t al. (20151. We use emcee 
(Foreman-lVlackey et al.|2013|, an implementation of the 
atfaneu nvariant ensemble sanipler of [Goodman fc Wearej 
(20101 to sample from the posterior probability distribu¬ 
tion tor the stellar—limb darkening coefficients, mass, ra¬ 
dius, and effective temperature—and planetary—radius, 
period, phase, impact parameter, eccentricity, and argu¬ 
ment of periapsis—parameters, conditioned on the light 
curve and t he measured stellar properties . 

Following Foreman-Mackey et al. (20151, the likelihood 
function that we use is marginalized over the weights of 
the “ei gen light curves” in the linear systematics model. 
Unlike Foreman-Mackey et al. (2015), we include an em¬ 
pirical Tjaussian'^rior'on~tIieweigEts determined by ro¬ 
bustly computing the distribution of weights across the 
full set of Campaign 1 light curves. This prior miti¬ 
gates the incorrect detection of false signals ind uced by 
stellar variability—as discussed below in Section [5.2) — so 
we exclude these candidates (EPIC 201929294 and EPIC 
201555883) from the tables of results. 

In this analysis, we assume the dilution caused by ad¬ 
ditional stars contributing flux into the aperture is neg¬ 
ligible for nearly all systems. Given the location of the 
Campaign 1 field at a high Galactic latitude, we expect 
low contamination by background giants. Nevertheless, 
this assumption may not be valid for all systems. Any 
contamination unaccounted for, as may happen if any of 
these stars are actually unresolved binaries, would cause 
us to underestimate the radii of any planets we detect. 
Therefore, high-contrast adaptive optics imaging of any 
systems should be obtained before these planets are used 
in population inference studies. The planet parameters 
measured by this analysis are listed in Table 3. 

4. FALSE POSITIVE ANALYSIS 

There are many scenarios which can cause an astro- 
physical false positive, where an EB star masquerades 
as a transiting planet. The most common scenarios 
are if (a) it is a highly grazing eclipse, or (b) the bi¬ 
nary system shares a photometric aperture with a sig¬ 
nificantly brighter star, resulting in a diluted eclipse 
depth. When possible, such astrophyscial false positive 
scenarios are traditionally ruled out by detailed follow¬ 
up observations, often a combination of high-resolution 


imaging and radial-velocity measurements. However, 
the Kepler mission, with its thousands of planet candi¬ 
dates around mostly faint stars, necessitated a paradigm 
shift—a move toward probabilistic interpretation of tran¬ 
sit signals, rather tha n comprehensive follow-u p of each 
individual candidate (Morton fc Johnson||2011 1. 


Morton (20121 presented an automated method to cal- 
ciilate the probability that a planet candidate might be 
caused by an astrophysical false positive. This method 
uses Galactic population simulations to determine the 
distributions of possible false positive scenarios, compar¬ 
ing the typical light curve shape of each to the data. It 
then combines this information with observationally mo¬ 
tivated prior assumptions about the populations of field 
stars, the properties of multiple star systems, and the 
occurrence rate of p lanets as determined from Kepler 
(Fressin et al. 2013|, in order to determine the proba- 


bility that the observed signal may be a false positive. 
Similar in spirit to other published methods of proba- 
bilistic validat ion, such as BLEN DER (Torres et al. 20111 
and FASTIS (Diaz et al. 2014), it has the advantage of 
being computationally less demanding and fully auto¬ 
mated, and thus easily applied in batch to a large number 

of candidates. __ 

In this work, we use vespep^ (Morton 2015bI, a new 
pu blicly available , general-purpose implementation of 
the Morton (20121 procedure, to calculate FPPs for each 
of these K2 candidates. The following constraints on 
false positive scenarios are imposed: 

• A chance-aligned EB system may reside anywhere 
inside or within one pixel of the photometric aper¬ 
ture of the target star. In creating a light curve for 
each star, we define photometric apertures rang¬ 
ing from 10 to 20 arcsec for each star, as defined 
in Table 4. Given the 6 arcsec point spread func¬ 
tion (PSF) of the Kepler telescope, we allow for 
the possibility that companions falling just outside 
of our aperture (within one pixel) may contribute 
to the light curve, possibly causing a false positive 
event. The search for such companions is discussed 
in ^ 


• The maximum allowed depth of a potential sec¬ 
ondary eclipse event is the most significantly de¬ 
tected signal at the same period of the planet can¬ 
didate, on ce th e primary transit is masked out (dis¬ 
cussed in 54.1). vespa does not allow for the pos¬ 


sibility of secondary eclipses larger than those ob¬ 
served in the K2 light curve for each star. 

• Blended stars must be allowed by the available 
adaptive optics an d ar chiv al im aging data (dis¬ 
cussed in detail in 54.2 and 4.3). vespa only con¬ 


siders stars below the detection threshold for the 
AO imaging, which is a position-dependent value 
following a calculated contrast curve for each star. 

Each of these scenarios is an astrophysical eclipse, 
caused by one object passing in front of another, blocking 
some fraction of the total light. The calculations here do 
not include the possibility that each signal is caused by 

http://github.com/timothydmorton/vespa 
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an instrumental artifact in the data or some other astro- 
physical event, such as stellar activity, masquerading as 
planet transits. 

Table 5 summarizes the results of these calculations, 
presenting the relative probability for each candidate 
to be caused by any of three false positive scenarios: 
an undiluted EB, a hierarchical triple eclipsing binary 
(HEB), and a chance-aligned background(/foreground) 
eclipsing binary (BEB). 

Six of the presented candidates have FPP >90%; these 
are considered to be likely false positives. On the other 
hand, 24 candidates have FPP < 1%. Three of the tran¬ 
sit signals might plausibly be caused by contamination 
by detected st ellar companions within the photometric 
apertures (see H.3l, so we keep these as candidates. 

This leaves 21 candidates that we statistically validate 
as planets, including fo ur that have been previously iden- 
tified in the literature ( Armstrong et al(]|2015 Crossfield 


et al. 20151. So in total, ot the 36 candidates, 21 are 


secure planets, 17 of which we validate here for the first 
time. 

We emphasize that the majority of these validations 
rely solely on the transit photometry and SDSS data, 
with follow-up imaging only obtained for seven of the 
31 targets. This demonstrates the utility of the vespa 
tool, which will be crucial to interpreting future candi¬ 
dates detected by K2, TESS, and PLATO and priori¬ 
tizing follow-up observing efforts. We show the transit 
signals in Figure 


4.1. Secondary Eclipse Observations 

One of the definitive signatures of a false positive bi¬ 
nary star system masquerading as a transiting planet is 
the presence of a secondary eclipse. While a nondetec¬ 
tion of a secondary does not exclude the possibility of a 
binary system (the orbit may be eccentric, or the com¬ 
panion too faint for a secondary eclipse to be detectable 
in the noise), such a nondectection reduces the probabil¬ 
ity of each of the EB false positive scenarios. 

To attempt to eliminate each EB scenario, we first 
search each K2 light curve to determine which secondary 
eclipse signals are not allowed by the data. We mask the 
transit signal of the planet in question and search for the 
most significant signal at the same period. Such a scheme 
does not assume circular orbits: we return the most sig¬ 
nificant signal at any phase, not only at the midpoint 
between consecutive transits. 

We report these maximum allowable secondary eclipse 
depths in Table 5. These values are used by vespa as 
limits on the allowable secondary eclipse. Any models 
that cause a larger event, such as a background EB con¬ 
sisting of two equal-mass stars in a circular orbit, can be 
excluded by the data. We note that with the exception 
of K2-19c, all systems with a maximum eclipse depth 
of at least one part per thousand have FPPs of 0.866 
or larger. The exception, K2-19, is a two-planet system 
with the two planets near a 3:2 period commensurability, 
so in this case the “secondary” is actually the transits of 
the other planet. 


4.2. Adaptive Optics Imaging 

We obtained high resolution images of seven stars 
with the Palomar High Angular Resolution Observer 


(PHARO) infrared detector (Hayward et al. 2 001|| be¬ 
hind the P ALM3000 adaptive optics system (Uekany 
et al.||2013l at the Palomar 5.1 m Hale telescope on the 
nipits'br2D15 February 3 and 4 UT. Sky conditions were 
mostly clear with light cirrus and «1'.'0-1'.'3 seeing on 
both nights. We used the smallest plate scale of 25 mas 
pixel”^ which resulted in a field of view of 25'.'6x25'.'6 
across the 1024^ pixeP array. All observations were ob¬ 
tained with the 32x pupil sampling mode, resulting in 
Strehl ratios of «20-30% in Ks for our P=ll-13 mag 
targets as measured by the Strehl monitor at the tele¬ 
scope in real time. We obtained unsaturated dithered 
frames of each target in ATs-band with typical integra¬ 
tion times of 2-10 s. Except for EPIC 201828749 and 
EPIC 201546283, which had nearby candidate binary 
companions, we also acquired deep saturated images (5- 
10 frames at 60 s each) to search for fainter companions. 

Images were r egistered and contra st curves were gen¬ 
erated following Bowler et al. (2015). For the saturated 
data, the star’s position in each image was found by 
masking the saturated region and fitting a 2D bivari¬ 
ate Gaussian to the PSF wings. Contrast curves for the 
median-combined image are calibrated using the unsat¬ 
urated frames. The typical sensitivity is 6.5-7.5 mag 
at 1". The images were astrometrically calibrated using 
dithered observations of the Trapezium cluster centered 
on 6^ Ori C taken on 2015 February 3 UT. Based on the 
reference astrometry for p a irs of stars in the field from 
McCaughrean & Stauffer (1994), we measure a plate 
scale of 25.2 ± 0.4 mas pixel”^ and north orientation of- 
0?2 ± 0?3. Since this latter value is consistent with being 
aligned with the detector columns, we adopt a value of 
0?0 ± 0?3 for this work. Relative photometry of nearby 
stars is carried out using aperture photometry with an 
aperture radius of 12 pixels (0'.'3). For EPIC 201828749, 
we also acquired J- and i7-band images. Astrometry and 
photometry is derived separately for each image, and the 
mean and standard deviation of these measurements is 
adopted for our final values listed in Table 4. 

Images for all systems AO data was obtained for is 
shown in Figure]^ while contrast curves showing the 5cr 
limits for detection as a function of orbital separation are 
given in Figure 

4.3. Known Background Stars 

The PHARO AO system has a field of view of 25 arcsec. 
Each K2 pixel is a square, 3798 on a side. A background 
EB within a few K2 pixels of our target stars could mimic 
a transit signal inside our aperture while evading detec¬ 
tion by PHARO. Such wide EBs should appear in seeing- 
limited ground-based surveys. 

To investigate the possibility that such wide compan¬ 
ions exist, we query the ninth data release of the Sloan 
Digital Sky Survey (SDSS DR9,|Ahn et al.||2012|). For 
each target, from the depth of tne observed transit we 
determine how bright a background object must be to 
cause the event if the background object were an equal 
mass totally eclipsing binary. We then search for all stars 
within 25" that are within this brightness limit relative 
to the candidate host star. All apertures we use in our 
K2 analysis are smaller than 20" so this search should 
encompass the region where possible background con¬ 
taminants could reside. Of the 31 stars in our sample, 
eleven have such a companion, plus one detected in AO 
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Fig. 2.— Phase-folded K2 photometry for all planet candidates analyzed in this paper. Each is the product of a fiducial noise model, in 
which the median systematic has been removed for illustrative purposes. The systems which we validate as transiting planets are labeled 
in blue. The systems which we confirm as false positive events are labeled in red. The systems which we leave as candidates are labeled in 
black. Red curves outline the median transit model for each candidate system. 


imaging. 

Unlike the original Kepler field, the field for K2 Cam¬ 
paign 1 is well out of the Galactic plane, so the rate of 
giant, distant background stars is significantly lower. We 
include all potential contaminants in Table 4. We vali¬ 
date or eliminate each of these as a possibility based on 
the transit shape. For example, the events near EPIC 
201546283 could only be caused by a background binary 
if the background object was a completely eclipsing sys¬ 
tem (so that the eclipse depth was 50%). In this case, 
the transit would be V-shaped. Since it is not, the back¬ 
ground object likely does not cause the transit event. 

In Table 4, the “maximum depth” column represents 
the maximum observed “transit” depth if the transit were 
actually caused by a total eclipse of the hypothetical 
background binary system, inducing a 50% flux decre¬ 
ment in the background star’s apparent brightness. 

The photometric apertures used to detect these can¬ 
didates range in radius from lOl'D to 19'.'9. In order to 
be a plausible contaminant, any companion star must 
be either within this aperture or just outside but bright 
enough for signficant flux to leak in. Evaluating each of 
the systems listed in Table 4, we judge that we cannot yet 
rule out contamination as a potential source of the transit 
signal for four candidates: 201295312.01, 201403446.01, 
201546283.01, and 201828749.01. Despite receiving low 
FPP scores from vespa, we list these systems as candi¬ 


dates in Table 5, rather than planets. Further updates to 
the vespa code will allow consideration of “specific” false 
positive scenarios; that is, scenarios that correspond to 
actually detected stars such as these, rather than hypo¬ 
thetical background or bound companions. 

The candidates with identified companions that we 
judge to not be plausible sources of potential contami¬ 
nation are the following: 

• K2-13b (201629650.01)— The companion to this 
star is 1773 from the EPIC target. As this is out¬ 
side the aperture (radius 1579) and the background 
star is not particularly bright, we rule out contam¬ 
ination for this system. 

• 201702477.01— The companion to this star is 
12715 from the EPIC target, and the aperture size 
is 1070. In addition, the maximum depth in this 
system is almost identical to the transit depth. For 
these two reasons we rule out contamination in this 
case. 

SDSS is 95% complete at r = 22.2 mag and the tele¬ 
scope has a PSF of 174. For the purposes of the vespa 
calculation, we thus treat nondetection in SDSS data as 
providing a contrast curve at wide separations down to 
a limiting magnitude of r = 22.2 mag. 

4.4. Archival Imaging 
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Fig. 3.— Adaptive optics images for the seven stars observed with high-contrast imaging. The main frame for each single system shows 
the deep, saturated image. The inset for each single system shows a shallower, unsaturated image to better identify companions at close 
projected orbital separations. For the two systems with imaged companions, EPIC 201546283 and EPIC 201828749, only unsaturated 
frames are collected. The pixel scale is 070252 per pixel. Each subplot is a square 400 pixels on a side and each inset is a square 100 pixels 
on a side. All subplots, including insets, are plotted on the same scale. 


For the stars with AO nondetections, there is still the 
possibility that a background binary could be positioned 
directly behind the target star, evading detection. The 
probability is small, given the CK'l diffraction limit of the 
Hale Telescope at 2 /rm, but nonzero. While the vespa 
calculations quantify this probability for this to occur, 
we can also rule out the possibility of such chance align¬ 
ments, down to a certain contrast, with archival imaging 
data. 

Five of the stars in our sample have proper motions 


larger than 50 mas yr—1, so they have moved across the 
sky by > 21'5 since they were imaged during the first 
Palomar Observatory Sky Survey (POSS) in the 1950s. 
To rule out background companions, we download data 
from the POSS I and II surveys, which imaged these 
targets in 1952-1955 and 1989-1998, respectively. We 
also download data from the Sloan Digital Sky Survey, 
which imaged these fields between 2000 and 2009. As 
shown in Figure we do not detect any background 
targets at the present-day location of any of these stars 























Fig. 4.— Scr contrast curves for all systems with AO nondetec¬ 
tions. For all systems, we can exclude the possibility that a com¬ 
panion at a given AKg exists. From our known transit depths, we 
can then rule out significant parameter space in which an eclipsing 
binary could reside and mimic a transit signal. 

in any of these images. 

For these targets, we can extend our contrast curves 
to zero present-day orbital separation and rule out the 
possibility that these transit events are caused by a back¬ 
ground EB. By combining present-day seeing-limited 
photometric survey data, adaptive optics imaging, and 
archival photometry, the only stellar companions we 
would not detect would be those that are gravitationally 
bound to the target star and positioned in their orbits so 
that their projected separation is smaller than the diffrac¬ 
tion limit of the Hale Telescope. Such an alignment 
would require the orbital inclination of the binary to be 
nearly 90° and the phase vu + 9 « 7r/2 or 37r/2. While 
we cannot fully rule out this possibility, the vespa cal¬ 
culations confirm that its probability is negligibly small. 


4.5. Tillinghast Reflector Echelle Spectrograph (TRES) 
Radial Velocities 

We observed K2-18 on 2015 February 04 and 25 UT 
with TRES on the 1.5 m Tillinghast Reflector at the Fred 
L. Whipple Observatory. These dates were chosen to be 
near the times of largest RV variations, corresponding 
to phases of 0.72 and 0.32 relative to the time of tran¬ 
sit. The spectra were taken with a resolving power of 
R = 44,000 and integration times ranging from 2800 to 
3600 s, resulting in S/N between 17 and 29 per resolution 
element. _ 

The spectr a were extracted as described in [Buchha^^ 
et al. (20101. The relative RVs were derived by cross- 
correlating the spectra against the strongest observed 
spectrum (in this case, the first) over the wavelength 
range 4700 - 6800 A. We selected 19 echelle orders in 
the analysis, being careful to reject orders with telluric 
absorption lines, fringing in the far red and those with 
very low SNR in the blue. 

The two observed spectra have RVs that differ by 47 ± 
42 m s“^. If the RVs were caused by a stellar companion, 
the RV shift between these observations would be on the 
order of km s“^. Therefore, we can rule out any stellar- 
mass companions that would be able to create this transit 
signal. 


5. POTENTIALLY INTERESTING SYSTEMS 
5.1. A Mini-Neptune with Earthlike Insolation 

The planet orbiting K2-18 may be an interesting target 
for atmospheric studies of transiting exoplanets. 

By combining archival and modern seeing-limited data 
with adaptive optics imaging, we can exclude the possi¬ 
bility these transit events are caused by a background 
EB. The apparent transits must be caused by an ob¬ 
ject co-moving with K2-18; radial velocities eliminate the 
possibility the companion is nonplanetary. Therefore, we 
confirm the planetary nature of this system. 

This star is an M2.8 dwarf at a distance of 34±4 pc. Of 
our pla net candidate hosts, o nly K2-3 (originally discov¬ 
ered by Crossfield et al. 20151 is brighter in RT- band. This 
star is only 0.1 m ag fainter in K than GJ 1214 (Charbon- 


neau et al. 20091. Due to the relative brightness of the 
host star, this target is likely to become a prime target 
for atmospheric characterization studies and is ideal as a 
target for future space-based missions such as JWST. 

The planet is slightly smaller than GJ 1214b, but un¬ 
like that planet, K2-18b is not highly irradiated. Instead, 
it is at a reduced semimajor axis a/i?*= 83.8 ± 9.0. Its 
equilibrium temperature is then, assuming zero albedo, 
Tf^q = 272 ± 15 K, meaning its bulk insolation is 94 ± 21 
percent that of the Eart h’s. Although the planet is likely 
too large to be rocky (|Rogers 20151, its atmosphere is 
likely to be the focus of many future observations, pro¬ 
viding a cool analogue to the highly irradiated planets of 
a similar size found by Kepler. 

5.2. Other Sources of False Positives 


The method of Foreman-Mackey et al. (20151 assumes 
that all variability in the light curves are caused by ei¬ 
ther the motion of K2, in which case the variability is 
shared by all stars, or transits of planets, in which case 
the variability is intrinsic to only one star. This assump¬ 
tion breaks down for extremely spotted stars where the 
astrophysical variability is larger than the instrumental 
magnitude. In that regime, the starspot modulations 
can be incorrectly fit by the systematic model, causing 
spurious transits to appear. This appears to be the case 
with EPIC 201929294, which has coherent starspots that 
appear to have the same rotation period as the transit 
period reported previously. Because the starspots are so 
periodic and coherent, these spurious transits were falsely 
identified as a planet candidate; we consider that system 
a false positive in this work. 

The candidate object possibly orbiting EPIC 
201555883 has a period, time of transit, and tran¬ 
sit duration consistent with EPIC 201569483. Such 
effects are n ot uncommon in Kepler data. |Coughlin] 
et al. (20141 identify 685 KOIs as false positives and 
outline four physical reasons why these anomalies may 
occur. EPIC 201555883 is a unique case in that it does 
not appear to fall under any of these cases. It falls 
on module 23, while EPIC 201569483 is on module 8, 
neither 180° away from nor on the same column as 
this candidate. Moreover, there is not any evidence of 
a mechanism that could cause a third star to induce 
both the appearance of a 7% eclipse on one module and 
an additional anomalous transit event on a different 
module. Instead, this candidate could be a false positive 
caused by a different systematic mechanism. 
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Fig. 5. — Archival maging for the five highest proper motion targets in our sample. In all cases, there are no background objects directly 
behind the present day location of the target (red circle) that could be missed by the AO observations. Modern SDSS imaging can also 
rule out wide companions that may have been missed at wide separations, beyond the AO field of view, such as the companion which can 
be seen in the images of K2-10. All figures are aligned such that north is up and east to the left. All subplots are on the same scale. 
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Foreman-Mackey et al. (20151 modeled the systematic 
enects in the K2 light curves using a linear combination 
of “eigen light curves” generated empirically by running a 
principle component analysis on the light curves of every 
star. This means that the training set includes the light 
curves for variable stars, EBs, and even transiting plan¬ 
ets. Again, this star has significant variability caused by 
starspots. In this case, the fitting procedure tries to ac¬ 
count for stellar variability using the eigen light curves. 
This overfit gives undue weight to eigen light curves that 
include the transits of EPIC 201569483, causing this spu¬ 
rious transit to occur. Again, we consider this system to 
be a false positive. As stated in Section 3, by including 
an empirical Gaussian prior on the weights for the eigen 
light curves in the linear systematics model, the signals 
observed on EPIC 201555883 and 201929294 are miti¬ 
gated, suggesting such a scheme should be employed in 
searching for planet candidates in future campaigns. 

The problem of over-fitting stellar variability using 
eigen light curves can also be solved by adding a stel¬ 
lar activity model to our fitting procedure. In this case, 
the spacecraft motion could be fit simultaneously with a 
model of starspot modulation, astroseismic oscillations, 
and planet transits. Such a model is currently under 
development (Angus et al. 2015, in preparation) 

5.3. Multiple Planet Systems 
Five of the systems reported by [Foreman-Mackey et al.] 


(20151 have more than one transiting candidate. One of 


by Crossfield et al. 

(2015 

(Armstrong et al. 2 

015|, 


Another of these is K2-19 


mensurability. The remaining three are all representa- 
tive of the multiple-planet systems observ ed by Kepler 


(Lissauer et al. 2011 Fabrycky et al.|[2014 l. Two of the 


systems are near a period commensurability and all three 
consist of mini-Neptune sized planets. 

We do not detect any significant transit timing varia¬ 
tions (TTVs) in any of these systems from the K2 data 
alone. K2-5 would be expected to have a TTV period 
of 117 days, but is likely too far from commensurability 
to have an observable TTV signal. K2-8 is expected to 
have a TTV period of 234 days, so this system may be a 
candidate for additional follow-up to constrain the sys¬ 
tem masses dynamically. The transiting planets orbiting 
K2-16 are near a 5:2 period commensurability. There 
is no evidence from Kepler of an abundance of planets 
near this period ratio, and so this may be coincidence. 
Follow-up observations may be warranted to search for 
an additional planet in this system forming a resonant 
chain, similar to those observed around other stars (e.g. 
Swift et al.||2013 Campante et al.||2015 |. 


5.4. Systems Orbiting Bright Stars 

One of the primary goals of K2 is the detection of tran¬ 
siting planets around bright stars that can be followed up 
from the groun d or with future spac e-based observatories 
such as JWST ( Howell et al.|[2014 l. Of our sa .mple. two 


systems orb it stars with K < \) mag: K2-3 (Crossfield 
et al. 2015|) and K2-18. An additional planet candidate 
may orbit EPIC 201828749, a star with K = 9.93 ± 0.03 
mag. These targets are ideal for ground-based followup 


and may be useful targets for Spitzer and JWST to probe 
planetary atmospheres. 

6. RESULTS AND DISCUSSION 

We have presented stellar pa rameters for all planet 
candid ates systems identified by [Foreman-Mackey et ah] 
(20151. We statistically validate 21 of the 36 candidates 
as bona fide planets, and we identify 6 as false positives, 
including two systematic false alarms. Of the planets, 4 
have been previously validated in other works, while 17 
are validated here for the first time. The systems not 
validated as planets or false positives remain as planet 
candidates. 

Enabling much of this analysis are two new open-source 
Python packages: isochrone^^ which we use to infer 
posteriors on physical stellar properties based on fitting 
theoretical stellar models to observed data; and vespcp] 
a new implementation of the Morton (2012 1 transit false 
positive analysis scheme. Both of these packages will 
continue to be useful in future analysis of transit candi¬ 
dates where comprehensive follow-up observations may 
be unavailable. 

The isochrones package uses the nested sampling 
scheme MULTINEST to capture the true multimodal na¬ 
ture of the posteriors. Using an MCMC algorithm in¬ 
stead can cause only one peak in the posterior distri¬ 
bution to be sampled. If the photometry is consistent 
with both a star on the main sequence and the sub¬ 
giant branch, an MCMC technique could cause one of 
these peaks (likely the subgiant possibility) to be missed, 
leading to an underestimation in the likelihood of sub¬ 
giant stars and and underestimation of the uncertainties 
of both the stellar and planetary parameters. 

With the exception of one object, all of the stellar pa¬ 
rameters are derived from comparing photometric obser¬ 
vations to the Dartmouth stellar evolution models. As a 
result, both the stellar and planet parameters are subject 
to systematic biases induced by discrepancies between 
the models and reality. 

The planets we confirm in this paper, like the planets 
found in the original Kepler mission, span a wide range 
of parameter space. They are at distances ranging from 
34 to 700 pc, have radii ranging from 1.3 to 5.3 i? 0 , 
and orbit with periods ranging from 5.0 to 50.3 days. 
Like the original mission, we find significantly more small 
planets than large planets, as expec ted from the radius 
distributions measured from Kepler ([How ard et al.|2012 
Fressin et al. pot's Morton fc Swift[[2014| . 


Unlike the origihal mission, however, we find that 
nearly all of our confirmed planets are around stars less 
massive than the Sun. This difference is a result of both 
the Campaign 1 field and the target selection process. 
Campaign 1 is at a significantly higher Galactic lati¬ 
tude than the original Kepler mission, meaning there is a 
much lower number density of targets at large distances. 
As massive stars at kiloparsec distances are relatively 
less likely to exist in Campaign 1 than near the Galac¬ 
tic plane, the pool of targets that could be selected for 
Campaign 1 contains a larger fraction of subsolar stars. 

Low-mass stars, particularly M dwarfs, are also a spe¬ 
cific focus of the K2 mission. One of the primary goals 

http: //github. com/timothydmorton/isochrones 
http://github.com/timothydmorton/vespa 
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of the Kepler mission was to “determine the abundance 
of terrestrial and larger planets in or near the habitable 
zone of a w ide variety of spectral types of stars” (Batalha 
et al.|2013 1 However, ~ 70% of Kepler’s target stars had 
masses within 20% of the Sun’s, while 70% of the stars 
i n the Galaxy hav e less than 50% the mass of the Sun 
(Brown et aljj2011 1. K2 will fulfill the promise of Kepler, 


with the goal ot providing a yield of small planets around 
bright, small stars t o facilitate follow-up measurements 
(Howell et al. 20141. This is clear from the K2 target 
selection process, with thousands of K and M dwarfs be¬ 
ing selected in each campaign. Based on these plans, we 
expect that K2 will detect hundreds of planets during 
its lifetime, with the majority being mini-Neptunes and 
super-Earths around stars less massive than the Sun. 
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TABLE 2 

Stellar Properties for all Objects of Interest 


EPIC 

Alternate 

RA (J2000) 

Dec (J2000) 

Mass 

Radius 

Teff 

[Fe/H] 

Distance 


Name 

(Degrees) 

(Degrees) 

(Mq) 

{R@) 

(K) 

(dex) 

(pc) 

201208431 

K2-4 

174.745639 

-3.905585 

0 63+0.03 
o'-oo_Q Q3 

U.DU_o 02 

419714 ! 

—0 12+0-10 
'^•-‘-^-0.12 

2181!! 

201257461 

— 

178.161110 

-3.094936 

-1 Kn+0.04 
-*-•0^-0.02 

iu.yD_o 93 

514ll!® 

— 0 21“*‘9-01 

16511!!! 

201295312 

— 

174.011629 

-2.520881 

1 07+®'®^ 
-‘-•O'l _0.07 

-‘-•O'^-O.il 

59891!°° 

—0 

'^•O'^-O.lS 

3311!! 

201338508 

K2-5 

169.303502 

-1.877976 

n ^^3+0.01 

0 52+0-01 

0-oz_o 01 

41021!° 

_Q t:-i+0.04 

^•O-L_0 06 

1811 ! 

201367065 

K2-3 

172.334949 

-1.454787 

n c;3+o o2 

0 52+0-02 
0-oz_o 02 

3951111 

—0 30+0-07 

42l° 

201384232 

K2-6 

178.192260 

-1.198477 

0 97+0 07 
-0.07 

o-yo_o 09 

58501!° 

— 0 14+0-17 

0'-l^-0.20 

343l°! 

201393098 

K2-7 

167.093771 

-1.065755 

0 97+0 06 
-0.06 

o-yo_o 08 

5772 I!! 

-0.071°;!° 

4331!° 

201403446 

— 

174.266345 

-0.907261 

-‘-•'^-‘--0.06 

1 1 9+0-26 
-‘-•-‘-^-0.14 

6445lf!4 

-0.501°;!° 

3621°° 

201445392 

K2-8 

169.793666 

-0.284375 

0 79+0 03 
'^-0.04 

0 74+0-02 
O-^^_0.03 

4890111 

1^-0'4_0.13 

4051!! 

201465501 

K2-9 

176.264467 

0.005301 

0 24+0 05 
'^•^^-0.03 

0 25+0-04 
O-^o_0 03 

3468+!! 

— 0 4b"*"9-19 

661!° 

201505350 

K2-19 

174.960319 

0.603575 

n ra+O-04 


55191!° 

—0 27+0-10 

29ll°° 

201546283 

— 

171.515164 

1.230738 

0.89tJ;J| 

0.881^1^ 

5422l!°34 

— 0 09+0-51 

2511°!°° 

201549860 

— 

170.103081 

1.285956 

n -TQ+O-OS 

'^•''^-0.03 

u.Dy_o 02 

45231!° 

0 05+0-15 
'0-Ob_o.i4 

2491! 

201555883 

— 

176.075940 

1.375947 

0 54+0-07 

0'-O^_0.01 

0 52+0-08 
0-oz_o 01 

44191!° 

-0.98l!|!! 

2891!° 

201565013 

— 

176.992193 

1.510249 

n 51 +0-15 

0'-01_Q Q3 

0 50+0-12 
d-OU_o 03 

39871!!° 

-0-441!!;!! 

5061!°° 

201569483 

— 

167.171300 

1.577513 

n cq+0.05 
^•O'^_0.05 

0 79+0-06 
o-^y_o.o5 

51921°° 

-0.09l°;!! 

1521!! 

201577035 

K2-10 

172.121957 

1.690636 

0 94+0-04 

'4-^^_0.06 

0 93+0.16 
u.yo_o 07 

56471°° 

O'-O'^-0.17 

2711!° 

201596316 

K2-11 

169.042002 

1.986840 

-I qr+0.04 

c 1 k+ 0.20 
0-lO_4 39 

5433l!°4 

— 0 12+0-01 

2019l!!28 

201613023 

K2-12 

173.192036 

2.244884 

1 01+®'®^ 

06 

1 01+0-27 
l-'Ol_0.09 

58001°° 

0.031°;!° 

294+!° 

201617985 

— 

179.491659 

2.321476 



37421°! 

-0.081°;!° 

lilt! 

201629650 

K2-13 

170.155529 

2.502696 


0.78l°;0®5 

56981!! 

—0 54+0-12 

290l!! 

201635569 

K2-14 

178.057026 

2.594245 

0 47+0-01 
O'-^^-O.Ol 

0 45+0-01 
0-^0_o 01 

37891!! 

r\ o't+0.03 
^•'J<_0.04 

2191® 

201649426 

— 

177.234262 

2.807619 

1 29+0-02 

Q 1 k+0.32 
®--‘-O_0.23 

50861!! 

—0 17+0-01 

25371°! 

201702477 

— 

175.240794 

3.681584 

n C7+0-06 
^•^^-0.06 

0 85+^-^^ 
d-oO_o 08 

561811° 

0'-^o_0.i8 

673l°! 

201736247 

K2-15 

178.110796 

4.254747 

n -To+O-OO 

O'- '^-0.03 

n rr+O- 06 
U-o8_0 03 

51311°° 

—0 4b"*"9-90 

4371!° 

201754305 

K2-16 

175.097258 

4.557340 

0'-O^_0.03 

O-O4i_0 03 

476 ll!! 

— 0 40"*"9-19 

3241!° 

201779067 

— 

168.542699 

4.988131 

0 91+0-03 

0 92+0-20 
^•y^_o.o7 

61661!! 

—0 54+0-07 

0'-O^-0.12 

188+f! 

201828749 

— 

175.654343 

5.894323 

0 74+1-06 

O'- '^-0.04 

0 71+9-04 

0-'1-0.06 

5552 I®! 

0'-oy_o.23 

1461!°°° 

201855371 

K2-17 

178.329776 

6.412261 

0 71+0-02 
0'-'J--0.05 

O-OO_0 03 

4320 I!! 

O 15+0-09 
'0-lO_Q 22 

1341° 

2019125521 

K2-18 

172.560461 

7.588391 

0 413+0-043 

u.^io_o 043 

0 394+0-038 

u.oy^_0 038 

35031°° 

u.uy_Q Q9 

34l! 

201929294 

— 

174.656968 

7.959611 

n -To+O.OO 
0'-^o_o 09 

Q yQ+0.04 

0-^u_o 08 

47861!° 

—0 16+0-22 
l^-lO_0.34 

1971!° 


Note. — These values and their uncertainties are derived from MULTINEST analysis and the numbers are computed as the 0.158, 
0.500, and 0.842 posterior sample quantiles. The coordinates are retrieved directly from the EPIC. These data are available in 
interactive form at https://filtergraph.com/k2_planets_montet. 


^ Parameters inferred from spectroscopic observations. 
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Planet Properties for all Objects of Interest 
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form at https://filtergraph.com/k2_planets_montet. 

^ Parameters inferred from spectroscopic observations. 

^ Declared a false positive due to noise modeling systematics (see §5.2) 



16 


TABLE 4 

Detected Companions to Candidate Host Stars 


Primary 

Aperture^ 

(arcsec) 

RA2 

(J2000) 

Dec^ 

(J2000) 

Detection^ 

Separation"^ 

(arcsec) 

ArS 

(mag) 

Max Depth® 
(ppt) 

Observed Depth^ 
(ppt) 

201208431 

15.9 

174.748988 

-3.902146 

SDSS 

17.25 ±0.15'" 

5.90 ± 0.12 

5.6 

1.20 

201257461 

19.9 

178.164376 

-3.093431 

SDSS 

12.91 ± O.IS*" 

5.04 ± 0.03 

4.8 

30.54 

201295312 

11.9 

174.010158 

-2.522528 

SDSS/AO 

8.12 ± 0.09*" 

7.10 ± 0.10 

0.8 

0.30 

201338508 

15.9 

169.308176 

-1.873647 

SDSS 

22.92 ±0.07*" 

4.35 ±0.03 

9.1 

1.07 

201367065 

19.9 







1.26 

201384232 

13.9 

178.195303 

-1.192501 

SDSS 

24.14 ± 0.06*" 

5.93 ±0.03 

2.1 

0.68 

201393098 

15.9 







0.53 

201403446 

15.9 

174.267663 

-0.909645 

SDSS 

9.78 ±0.14*" 

4.56 ±0.08 

7.5 

0.23 

201445392 

13.9 







0.78 

201465501 

11.9 







2.83 

201505350 

19.9 







2.64 

201546283 

17.9 

171.515265 

1.229950 

SDSS/AO 

2.98 ±0.05"" 

5.87 ± 0.06 

2.3 

2.33 

201549860 

13.9 

170.097556 

1.288007 

SDSS 

21.21 ± 0.05*" 

2.26 ±0.03 

62.3 

0.80 

201555883 

10.0 







3.50 

201565013 

10.0 







45.8 

201569483 

19.9 







160 

201577035 

19.9 

172.118116 

1.687798 

SDSS 

17.19 ±0.12*" 

5.40 ± 0.03 

3.5 

1.44 

201596316 

15.9 







0.70 

201613023 

19.9 







0.42 

201617985 

15.9 







1.10 

201629650 

15.9 

170.158905 

2.502107 

SDSS 

12.30 ± 0.14*" 

5.98 ± 0.06 

2.0 

0.58 

201635569 

11.9 







9.43 

201649426 

19.9 







216 

201702477 

10.0 

175.238916 

3.678764 

SDSS 

12.15 ±0.12*" 

4.65 ±0.09 

6.9 

6.70 

201736247 

13.9 







1.21 

201754305 

11.9 







0.80 

201779067 

19.9 







84.9 

201828749 

11.9 

175.645724 

5.894714 

AO 

2.46 ± 0.04"" 

2.0 ±0.1"" 

137 

0.76 

201855371 

19.9 







0.99 

201912552 

13.9 







2.85 

201929294 

19.9 







13.56 


^ Defined aperture used to create the K2 stellar light curve. 

^ Position of imaged companion. 

^ Dataset used to detect the imaged companion. 

Distance between the primary K2 target star and companion, in the dataset in which the companion is detected. 

^ Difference in r-band magnitude between the primary K2 target star and the companion. 

® Observed “transit” depth if the imaged companion’s flux were fully contained in the aperture and if i t we re an equal-mass eclipsing binary, leading 
to an eclipse depth of 50%. This is the maximum possible false positive eclipse depth, as described in |4.3| 

^ Observed transit depth in the K2 dataset. If larger than the "max depth," this transit event cannot be caused by eclipses of the background star. 
^ Separation from AO imaging 
^ Separation from SDSS photometry 
^ Ar inferred from JHK relative photometry. 
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TABLE 5 

False Positive Probability Calculation Results 


Candidate 

^sec,max [ppt] 

ao?2 

PrEB 

PfBEB 

PfHEB 

/p" 

FPP 

Disposition 

201208431.01/K2-4b 

0.51 

_ 

< 10-* 

8.1 X 10“* 

< 10-* 

0.21 

8.1 X 10-* 

Planet 

201257461.01 

0.59 

- 

0.998 

1.7 X 10-3 

< 10-* 

0.00 

1.000 

FP 

201295312.01 

0.04 

Y 

1.4 X 10"* 

< 10-* 

< 10-* 

0.17 

1.4 X 10-* 

Candidate^ 

201338508.01/K2-5C 

0.63 

- 

< 10-* 

2.9 X 10-3 

< 10-* 

0.22 

2.9 X 10-3 

Planet 

201338508.02/K2-5b 

0.33 

- 

< 10“* 

1.7 X 10-* 

< 10-* 

0.22 

1.7 X 10-* 

Planet 

201367065.01/K2-3b 

0.15 

- 

< 10“* 

1.1 X 10-* 

< 10-* 

0.22 

1.1 X 10-* 

Planet'^ 

201367065.02/K2-3C 

0.67 

- 

< 10“* 

< 10-* 

< 10-* 

0.16 

< 10-* 

Planef^ 

201384232.01/K2-6b 

0.44 

- 

8.4 X 10-3 

< 10-* 

< 10-* 

0.07 

8.5 X 10-3 

Planet 

201393098.01/K2-7b 

0.52 

- 

< 10-* 

1.1 X 10-3 

< 10-* 

0.05 

1.1 X 10-3 

Planet 

201403446.01 

0.18 

Y 

4.8 X 10-* 

< 10-* 

< 10-* 

0.19 

4.9 X 10-* 

Candidate^ 

201445392.01/K2-8b 

0.26 

- 

< 10-* 

2.1 X 10-3 

< 10-* 

0.18 

2.1 X 10-3 

Planet 

201445392.02 

0.18 

- 

< 10-* 

0.019 

< 10-* 

0.21 

0.019 

Candidate 

201465501.01/K2-9b 

0.68 

- 

< 10-* 

5.8 X 10-3 

< 10-* 

0.21 

5.8 X 10-3 

Planet 

201505350.01/K2-19C 

2.69 

- 

< 10-* 

5.6 X 10-3 

< 10-* 

0.04 

5.6 X 10-3 

Planet'^ 

201505350.02/K2-19b 

0.70 

- 

< 10“* 

1.6 X 10-* 

< 10-* 

0.07 

1.7 X 10-* 

Planet'^ 

201546283.01 

0.15 

- 

7.0 X 10-* 

2.6 X 10-* 

< 10-* 

0.00 

9.6 X 10-* 

Candidate^ 

201549860.01 

0.18 

- 

< 10“* 

0.026 

< 10-* 

0.04 

0.026 

Candidate 

201555883.01 

0.94 

- 

- 

- 

- 

- 

- 

Fpb 

201565013.01 

1.69 

- 

0.783 

7.3 X 10-3 

0.063 

0.07 

0.853 

Candidate 

201569483.01 

2.06 

- 

0.822 

< 10-* 

0.174 

0.00 

0.996 

FP 

201577035.01/K2-10b 

0.14 

Y 

4.4 X 10-* 

< 10-* 

< 10-* 

0.07 

4.4 X 10-* 

Planet 

201596316.01/K2-llb 

0.45 

- 

< 10-* 

1.2 X 10-3 

< 10-* 

0.06 

1.2 X 10-3 

Planet 

201613023.01/K2-12b 

0.08 

Y 

< 10-* 

< 10-* 

< 10-* 

0.18 

< 10-* 

Planet 

201617985.01 

0.27 

- 

< 10-* 

0.012 

< 10-* 

0.18 

0.012 

Candidate 

201629650.01/K2-13b 

0.43 

- 

5.9 X 10-* 

2.0 X 10-* 

< 10-* 

0.13 

7.8 X 10-* 

Planet 

201635569.01/K2-14b 

0.79 

- 

< 10-* 

4.9 X 10-3 

< 10-* 

0.05 

4.9 X 10-3 

Planet 

201649426.01 

3.10 

- 

0.896 

< 10-* 

0.104 

0.00 

1.000 

FP 

201702477.01 

0.70 

- 

0.137 

1.2 X 10-3 

6.6 X 10-3 

0.05 

0.145 

Candidate 

201736247.01/K2-15b 

0.42 

- 

4.8 X 10“* 

2.1 X 10-* 

< 10-* 

0.19 

6.9 X 10-* 

Planet 

201754305.01/K2-16C 

0.65 

- 

1.0 X 10“* 

1.4 X 10-3 

< 10-* 

0.21 

1.5 X 10-3 

Planet 

201754305.02/K2-16b 

0.38 

- 

2.3 X 10-* 

9.9 X 10-* 

< 10-* 

0.19 

1.2 X 10-3 

Planet 

201779067.01 

1.97 

- 

0.968 

1.3 X 10-3 

7.2 X 10-3 

0.00 

0.976 

FP 

201828749.01 

0.39 

Y 

0.644 

3.8 X 10-* 

< 10-* 

0.01 

0.645 

Candidate 

201855371.01/K2-17b 

0.62 

- 

< 10-* 

8.7 X 10-3 

< 10-* 

0.01 

8.7 X 10-3 

Planet 

201912552.01/K2-18b 

0.47 

Y 

< 10-* 

< 10-* 

< 10-* 

0.21 

< 10-* 

Planet 

201929294.01 

3.12 

- 

- 

- 

- 

- 

- 

Fpb 


Note. — Results of the vespa astrophysical false positive probability calculations for all candidates. Likely false positives (FPP 
> 0.9, or otherwise designated) are marked in red. Candidates are declared to be validated planets if FPP < 0.01. EB, BEB, 
and HEB refer to the three considered astrophysical false positive scenarios, and the relative probability of each is listed in the 
appropriate column. Planets previously identified in the literature are marked. 

^ Maximum depth of potential secondary eclipse signal. 

^ Whether adaptive optics observation is presented in this paper. 

^ Integrated planet occurrence rate assumed between 0.7x and 1.3x the candidate radius 

^ Despite low FPP, returned to candidate status out of abundance of caution due to secondary star detection within or near 

photometric aperture. _ 

° Declared a false positive due to noise modeling systematics (see g5.2|l. 

^ Identified as planets by 
^ Identified as planets by 


Crossfield et al.| 

2015L 

Armstrong et al. 

(|2015ll 











